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ABSTRACT 

The modulation of galactic protons and He nuclei during 
the last solar cycle is analyzed according to Parker's theory. 

The mechanism of modulation remains essentially the same 
during several years of low solar activity (1961-1965). The 
modulation near solar maximum (1959) implies that the scale 
sizes of the magnetic inhomogeneities in the solar wind are 
reduced below the values at solar minimum. An adequate de- 
scription at solar maximum would require further refinements 
of the theory. The proton to He nucleus ratio outside the solar 
system is shown to be consistent with a value ~ 6, in a klnetic- 
energy/nucleon representation for the interval 50 to 1000 Mev/nucleon. 
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INTRODUCTION 

Since the discovery by Forbush [1954] of the inverse correlation 
of cosmic ray intensity with solar activity, many experimental and 
theoretical studies of the long term solar modulation effects on galactic 
cosmic radiation have been attempted* The present status of the subject 
has been presented in several recent reviews, [ Webber , 1964; Quenby, 
1964] as well as in some recent papers [Fichtel et al ., 1965; Freier 
and Waddington, 1965; Nagashima et al., 1966]* In this paper we analyze 
the available energy spectrum data obtained during the last solar cycle 
to see how far they conform to the predictions of the theory of the cosmic 
ray modulation caused by the solar wind [Parker 1963]. It is considered 
pertinent to do this as recent satellite measurements [Bonnetti et al ., 

1963; Ness et al ., 1964, 1966] have confirmed the existence of the solar 
wind, the solar origin of the interplanetary magnetic field and the con- 
vection of magnetic inhomo gene iues of the right size to modulate galactic 
cosmic rays. 

PARKER'S THEORY OF THE 11 -YEAR SOLAR MODULATION 

Parker [1963] has developed a theory of the 11 -year cosmic-ray 
variation in which the cosmic ray density at the earth's orbit is deter- 
mined by the steady state equilibrium between the inward diffusion of 
galactic cosmic rays and their outward convection by the solar wind. 
Assuming spherical symmetry Parker arrives at the following equation re- 
lating the cosmic ray density at the orbit of the earth fi e (/3) to the cosmic 

ray density outside the solar influence jjl (£): 

00 ' 
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^.(0) = /*»(£) eX P 



V(r) 
D(r, p) 


dr 


(i) 


where V(r) is the solar wind velocity and D(r, J3) is the diffusion coef- 
ficient, and the distance of the earth from the sun. In general 
D (r ,jS ) is a function of the radial distance r and the particle velocity 
/3c given by 

D(r, p) - /3c\(r, p) (2) 


where A.(r, p) is the mean free path for the diffusion process. For an 
arbitrary distribution of the scale size characterizing the inhomogeneities 
in the magnetic field 


1 

^(r, P) 



<?(%, p) p(l, r) d't 


(3) 


where o-(^, p) is the effective cross-section for the scattering of a 
particle with velocity /Sc by an inhomogeneity in the magnetic field with 
a scale size and p(^ f r) d'i is the spatial density of scattering centers 
with scale sizes between % and 'i + d'fc, at a distance r. 

We assume that r) is a separable function of 'i and r; 

r) = n(r) N(^) (4) 

where 

f N(^) d-t = 1 


0 
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We then have 


M(/3) = 




exp 



v. j 

In order to compare this modulation function with the observations, 
we need to specify a scattering cross-section aft, /3) and a distribution 
N(^)of the scale size of the inhomogeneities. 


Parker [1963] adopted 


2 for R < /£ 


aft,/3) = < 


ftVR 2 ) for R > l 


(6a) 


(6b) 


no£)= see-4 0 > 


(7) 


where R is the radius of curvature of a particle in a magnetic field B 
assumed independent of r: 


P - - A m ^2 fi . 1 

R "ZeB _ Z 0 (1-/3 2 ) 1/2 e B 


(8) 


where A and Z are the particle mass and charge numbers respectively 
and m 0 is the proton rest mass. 


With this choice of (3) and N(^) we can write 


- I for 




M(/3 ) = exp 


K 

/3 


(9) 
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M(/3) = exp 





for fi > fi 0 


where 


K 




3V(r) n(r) 


dr 


( 10 ) 


(ID 


R 


o ~ 


m o c2 

eB 


( 12 ) 


and 


^0 


(^ 0 Z/R 0 A) 

r i 1/2 

1 + (^ 0 z/r 0 A) 2 _ 


(13) 


In the appendix we write the modulation function M ( /S ) for the more 
general case where the magnetic field B is a function of the radial dis- 
tance r and N ('i ) is a general distribution function. The formula so 
obtained is presented for the limits 0 and fi -* 1, where the influence 
of the radial dependence of B and of the choice of the distribution N ) 
is made more explicit. 

EXPERIMENTAL DATA DURING THE PERIOD OF 
LOW SOLAR ACTIVITY 

Figure 1 shows the Mt. Washington neutron monitor counting rate 
(courtesy of Dr. J. Lockwood, University of New Hampshire) and the 
sunspot number during the last solar cycle. The arrows indicate the times 
of balloon flights from which the balloon data used in this paper were 
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obtained. The hatched region shows the time of the coverage provided by 
the IMP-I, II and III, and OGO-I satellites. 

Figure 2 shows the proton and helium nucleus energy spectra in a 
kinetic energy per nucleon representation. The helium nucleus flux was 
multiplied by 5. The curve labeled 1963 has the analytical form 

10 8 x * 8 

N(E) dE = 1U — dE (14) 

(E + 500) 4 

when E is in Mev. Equation 14 was constructed to fit the experimental 
data shown, and to have the asymptotic form N(E) aE" 2-5 at high energies 
[ Balasubrahmanyan et al., 1965 b ]. From this curve labeled 1963 
we constructed the curve lableled 1965 by assuming that AK 1965 = K J9l>5 y 19Ss 
= -0.2. In a like manner we constructed the curve labeled 1961, 

for AK 1961 = +0 * 4 - 

These data suggest that in the energy range of the observations 
and at a time near minimum solar activity, the solar modulation is 
predominantly velocity dependent. Studying the time variation of the 
low energy He nucleus spectrum Fan et al., [1965a ] , and 
Gloeckler [1965], arrived at a similar conclusion, for energies up to 
500 Mev/nucleon. The velocity ft q9 at which the modulation function 
changes from pure velocity dependence to a mixture of velocity and 
rigidity dependence, is a function (equation 13) of (4 q /R 0 ). The experi- 
mental data indicate that up to approximately 1 Bev/nucleon the modulation 
is only velocity dependent. From this we conclude that from 1961 
through 1965 the scale size of the inhomogeneities was such that 
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(* 0 /R 0 ) ^ 2 * For /? = 5y = 5 X 10‘ 5 gauss, this conclusion leads to 
l Q > 10‘ 2 AU. 

In addition to the low energy detectors, the satellites IMP-I, IMP-II, 
IMP-III and OGO-I carried GM counter telescopes of identical construc- 
tion. \Vhen the telescopes were calibrated with the sea-level cosmic- 
ray muon flux the counting rates agreed within 1%. Details of these GM 
counter telescopes are given elesewhere [Balasubrahmanyan et al., 
1965a ] . There is no detectable instrumental drift among the detectors 
from these different satellites and the maximum difference in the 
absolute counting rate was less than 5%. This was easily corrected by 
using data from periods when two or more satellites overlapped in 
time. Figure 3 shows the cosmic ray omnidirectional intensity above 
50 Mev measured by these GM counters from 1963 to 1965. 


Using the differential energy spectrum measured in 1963, we have 
calculated for each AK the expected ratio I 1955 /I 1963 of the integral 
intensity above 50 Mev in 1965 to the corresponding intensity in 1963. 
This calculated ratio is shown as a function of AK in Figure 4. The value 
of this ratio obtained from the GM telescopes falls within the hatched 
horizontal band shown in the figure, and therefore we conclude that the 
value of AK should be roughly between 0.18 and 0.Z0. This is consistent 
with the value 0.20 used to characterize the change in the low energy 
differential spectrum between 1963 and 1965. 
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EXPERIMENTAL DATA DURING A PERIOD OF HIGH 
SOLAR ACTIVITY 

During the period of high solar activity, the only suitable proton 
and He nucleus differential energy spectrum measurements were made 
by Webber and McDonald [1964], and Freier and Waddington [1965 ]. 

The energy spectra of protons and He nuclei (multiplied by 5) are shown 
in Figure 5. Below about 1 bev/nucleon the He nucleus spectrum is 
distinct from the proton spectrum. This situation is different from that 
at solar minimum where the spectra are essentially similar in a kinetic 
energy per nucleon representation. 

From Parker l s theory of the 11 -year solar modulation, with the 
choice of cr('fc) and N('C) given by equations 6 and 7, we have seen that 
the modulation is described by two different functions (equations 9 and 
10). Although equation 9 predicts a modulation which is the same for 
protons and He nuclei, equation 10 depicts a different modulation for 
these two components. Written in terms of the velocity as it is done in 
(10), the difference between the proton and He nucleus modulation lies 
in the factor (Z/A) 2 in the exponent. If the value of (^q/Rq) is sufficiently 
low (such that the transition energy corresponding to the velocity /3 Q is 
of the order of 100 Mev/nucleon), then we could expect that above this 
energy the He nuclei will be modulated less than the protons. Therefore 
the He nucleus intensity (after multiplication by 5) would be above the 
proton intensity and in qualitative agreement with the experimental 
observation. From the condition that fi 0 < 0.45 for He nuclei (correspond- 
ing to a kinetic energy /nucleon of 100 Mev) we found that in 1959 (^ 0 /R 0 ) < 1. 
For B = 5y this implies - 5 X 10 3 AU. 
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Although this conclusion about the value of Q can be reached by 
considering only the split between the proton and He nucleus spectra, a 
detailed fit of these spectra for various value of K and ^ 0 /R q is not 
possible. In order to investigate whether the origin of this difficulty 
lies in the assumption that the magnetic inhomogeneities have a scale 
size characterized by a single value 0 (equation 7), we have repeated 
the calculation using the following a and N (•£ ) : 


= r 

1 + (R/l) 2 


(15) 


N(4) = 1/(2t7) 1/2 -tjN exp 


(*-*o ) 2 

2^2 


(16) 


Equation 15 provides a continuous analytical expression for cr and 

in the limits of R « % and R » 't it tends to the values given by equa- 
tions 6a and 6b, Equation 16 represents a normalized gaussian distri- 
bution for N ("£) with a most probable value *t Q and a variance The 

normalization of N (Z) requires that 



With this choice of o-(^, J3) and N (Jt) we arrive at essentially the same 
results and conclusions reached from the simpler choice represented by 
(6) and (7). This is true for both periods of low and high solar activity. 
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PROTON TO HE NUCLEUS RATIO 

We have used the forms of cr(^, p) and N(/&) depicted by equations 
15 and 16 to investigate the proton to He nucleus ratio as a function of 
energy throughout the solar cycle. With this choice, the modulation 
function can be written as 

M(/3) = exp g (/?, a, <f, (18) 

where 


K' 


(277) 1 / 2 Nj V 


(19) 


g (A a, A V) = 


r i 


x 2 + (a<fy/?) : 


! ] 1 exp £- 


(x-1) 2 

2v 2 


dx (20) 


a =y, . V = ~j~~ and y = (1 -/3 2 )" 1/2 

Z 4 o % 

The proton to He nucleus ratio R pa (r e ) at the earth's orbit for the 
same kinetic energy /nucleon can then be written as a function of the 
velocity of the particle: 

V (0= R pa (°°) exp [-* (gp-ga)! (21) 

where 

gp = g(A 1, A T/), ga = g(A 2 , T)) 


and R^^oo) is the corresponding ratio outside the solar system* 
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The expressions g p and g a were integrated numerically from 
10 Mev/nucleon to 10 bev/nucleon for different values of the parameter 
£ . The parameter 77 was taken as 0.67. 

Figure 6 presents the ratio [R^ (r e )/R pa (ao)] as a function of the 
kinetic energy/nucleon for different values of K L and £~ l (= ^ 0 /R 0 ). 

K l is defined by 


K L = K' g (/3 = 0) (22) 

With this definition 

M(/3) - exp (~K l //3) (23) 

0-0 

for any 

The measured values of R p0 (r e ) for 1961, 1963 and 1965 yield 
values for (R pa (r e ) / R pa (oo)) all within the hatched area between curves 
A and C (Figure 6 ) providing we chose R pa (°°) = 6 . Figure 6 shows that, 
for the choice R ( r ^) = 6 , all the 1959 values for R pa (r e )/R pa (co) fall 
well outside the hatched area. 

We conclude that the data on the proton to He nucleus ratio from 
1959 through 1965, in the kinetic energy /nucleon interval of ~ 50 Mev 
to ~1 bev, are consistent with a constant proton to He nucleus ratio of 
6 outside the solar system. From Figure 6 we also note that the 1961-1965 
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data are consistent with (l 0 / R 0 ) > 2, and the 1959 with {l 0 /R q ) < 1. 
Although the proton to He nucleus ratio during 1959 can be fitted with 
several different combinations of K L and (^ 0 /R 0 ) (see Figure 6), the 
actual proton spectrum generated with any of these combinations fails 
to agree with the experimental spectrum. In this connection it is of 
interest to note that energy losses in the expanding solar wind (neglected 
in the treatment presented in this paper) may be important [Parker , 1965]. 
For large values of K, characteristic of solar maximum, the actual in- 
tensity decrease could be much larger than that predicted by the simple 
theory. As the experimental data are too scarce, it has not been pos- 
sible to draw any reliable conclusion regarding the role of energy losses 
in the modulation mechanism during the high solar activity period. 

PROTON AND HE NUCLEUS RIGIDITY SPECTRA 

Figures 7a, 7b, 7c and 7d show the proton and He nucleus (multiplied 
by 7.2) rigidity spectra during 1965, 1963, 1961 and 1959, respectively. 
During 1959 the proton and He nucleus spectra may be represented by 
the same curve, whereas during 1961, 1963 and 1965 there is a con- 
siderable split between these two components at low rigidities. This 
shows that during 1959 the modulation was essentially rigidity dependent 
down to about 600 Mv, indicating a small value of (^ 0 /R q ) near solar 
maximum. On the other hand, during 1961, 1963 and 1965 any rigidity 
dependent modulation could occur only for rigidities > 2 Bv, which 
indicates the large value of (-£ 0 /R 0 ) near solar minimum. 

McCracken and Rao [1966], analyzing the average cosmic-ray 
solar diurnal anistropy obtained by means of neutron monitor data, 
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arrive at the conclusion that the size and frequency of occurrence of 
the small scale irregularities in the interplanetary field did not change 
appreciably during the last solar cycle. We feel that this conclusion is 
not necessarily in contradiction with the views presented in this paper 
since the primary energies involved are different. In addition, the 
balloon and satellite data on which our conclusions are based refer to 
differential energy measurements, whereas the neutron monitor data 
reflect the integrated effect on particles with energy above ~1 bev. 

Since (^ 0 /R 0 ) changed by a factor of more than Z from 1959 to the 
years near solar minimum (1961-1965), it is important to relate this 
cosmic ray observation to direct solar observations. In this connection 
we note that Babcock [19593 has observed that, during the middle of 
1957, the polarity of the magnetic field near the south heliographic 
pole was reversed and that the field near the north pole was not seen 
to reverse until the end of 1958. There was a period of over a year when 
both the poles had the same polarity. Also, there was a six month period 
in 1957 when there was zero effective field. This unstable situation 
probably indicates a high overall turbulence giving rise to magnetic 
inhomogeneities with a wide range of scale sizes, including small values. 

The evidence presented by Hewish [1958] on the basis of scattering 
of radio waves from the Crab Nebula when it is occulted by the solar 
corona also tends to show that the general scattering increases near the 
solar maximum. Though these inhomogeneities may be on a different 
scale, it is likely that the spectrum of turbulence tends to include a large 
number of smaller scale sizes. 



* > 
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CONCLUSIONS 

1. During the period when the solar activity is low, in the interval 
20 Mev/nucleon to 1 bev/nucleon, protons and He nuclei are modulated 
in a velocity dependent manner in reasonable accord with Parker's 
theory. The intensities during 1965 and 1961 are related to the intensity 
during 1963by I 196S = I 19M X exp (+0.2/0) and I, 961 = I, 963 X exp ( -0.4/yS). 

2. The parameter (^ 0 /R 0 ) characterizing the magnetic turbulance 
during the solar minimum years 1961 - 1965 is equal to or greater than 2. 
During the period of high solar activity (1959) the existing data suggest 
that ('t /R ) is less than 1. 

3. The proton and He nucleus energy spectra near solar maximum 
derived from Parker's solar wind modulation theory, with a reasonable 
choice of parameters, cannot be made to fit the experimental observations. 
This suggests that energy losses in the expanding solar wind may be 
important near solar maximum. 

4. By an examination of the data in the interval 20 to 1000 Mev/nucleon 
for the years 1961, 1963 and 1965, it is concluded that the proton to He 
nucleus intensity ratio R (o°) outside the solar system is approxi- 
mately 6. This value is in close agreement with the universal abundance 
of these two elements. If this ratio increases at higher energies, as 

some measurements seem to indicate [Webber 1964], then there is a 


possibility that the proton and He nucleus spectra at the source are no 
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longer similar at higher energies. It is likely that a study of this ratio 
as a function of energy could give information on the mechanism of 
escape of these particles from the source regions. 
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APPENDIX 


In deriving Formula 5 of the text we have assumed that the magnetic 
field B could be considered constant, equal to some effective value B eff . 
Let us see now how this formula changes if we relax this restriction and 
let B be a function of the radial distance B(r). Then, for any distribution 
of the scale size of the inhomogeneities N(-£) we have for the modulation 
function 


M(0) = exp^-i 




''l 


cr(i , A r) N(^) d-C] dr 


(Al) 


where the cross section cr (■£ , fi, r) now depends also on the radial dis- 
tance r through its dependence on the magnetic field B(r). 


In order to see in more detail what influence B(r) and N(^) will 
have on M03) let us assume that 


a (-C, p, r) = -C 2 when P -> 0 


(A2) 


and 


a (t, P, r) = (^ 4 /R 2 ) when P - 1 


(A3) 


in agreement with (6a) and (6b). Then we obtain 


/•CD /* 00 

I cr(4, p, r) N(-£) d't = ■i 2 N(-t) d-t^ < ■i 2 > when P - 0 (A4) 

^ A ^ A 
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a(-i, /3, r) N('t) d't 



• l 4 N(-t) d'C ~~~J < > when /3~* 1 


( A5) 


and thus 


M(/S) - exp 
P-+ o 


/8 



3V(Q n(r) 
c 



M(/3) = exp 
/ 3 — 0 


(l-l 2 ) . 

/3 3 


z \ 2 <^ 4 > 

a) -5T 



3V(r) 

c 



n(r) dr 


( A6) 

(A7) 


where now R n = mc 2 /eB and B = B(r ). 

U e e x e / 


From this, we can see the effect that N(/£) and B(r) have on the 
low- and high-energy approximation formulas for M(/3) (equations A6 
and A7). The effect of choosing different forms for the distribution 
function N(^) is felt through the different values that < 'i 2 > and < 4 4 > 
will have in the low- and high-energy approximations. The choice of the 
function B(r) will not affect (A6) but introduce a variation in (A7) 
through the term 


f^-[W "<'><*' 


where different values for B(r) will result in different values for this 
integral. 
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If we now specify 


N( 4 ) 


(277) 1 / 2 l x Nj 


exp 


2^2 


then 


M(/S) = exp 
P— o 


*5 (1 + ^ + bij) f" lY( r )"< r ) dr 


exp 


_ K(1 +r? 2 + bT7) 

fi 


] 


and 


M(/3) = exp 

/ 3 - 1 


_ a-£) . ( zY . /M 


/ 3 3 


A/ \ R, 


^0 j^l+61? 2 +37j 4 +b (r)+ 57j 3 )J 


■f^[¥F 


n(r)dr 


= exp 


La- . (z\ 2 ( l o \ 2 


, A j \ R. 


Kj + 6t] 2 + 3i7 4 + b(r) + 5t7 3 )J 


where 


b = 


(277) 1 / 2 N exp (1/2 v 2 ) 


V =■ 


= -eg J” milDill dr , and K, = < 3 J" 3v <0 • p<Oj n(r) dr 

r r e 


(AS) 


<A9) 
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If in addition B(r) = constant = B e , then K x = K and equations A8 
and A 9 would agree with the low- and high-energy limits of (19) given 
in the text. 

On the other hand if 

NCC) = S(*-4 0 ) 


then 


M(/3) = exp 

y8-0 



M(/3) = exp 



(A10) 


(All) 


And again if B(r) = constant = B e , then K = K x and equations A10 


and All agree with (9) and (10) of this paper. 
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Fig. 1 


Fig. 2 


Fig. 3 


Fig. 4 


FIGURE CAPTIONS 

> Sunspot number and Mt. Washington neutron monitor rate during 
solar cycle 19. The arrows indicate the times of balloon flights 
from which the data used in this paper were obtained. The 
hatched region shows the time of satellite coverage. 

. Proton and He nucleus (x5) energy spectra measured by balloons 
and satellites during 1961, 1963, and 1965. The 1961 data are 
from Fichtel et al. [1965], Bryant et al. [1962], and Meyer and 
Vogt [1963]. The 1963 data are from the following sources: 
McDonald and Ludwig [1964], Fan et al. [1965a], and 
Balasubrahmanyan and McDonald [1964]. The 1965 data are 
from Balasubrahmanyan et al. [1965cJ. The lowest energy He 
nucleus data point exhibited for 1965 might have a contribution 
associated with long-lived solar streams [Fan et al, 1965b]. 

. The total integral intensity of cosmic rays above 50 Mev as 

measured by GM counters on IMP-I, II, III and OGO-I satellites. 

. The ratio of the 1965 to the 1963 total cosmic ray 

integral intensity above 50 Mev as a function of the increment 
AK in the parameter K of the solar wind modulation theory. 

The experimentally observed ratio Ii 9 65 /li 963 with! 11 the 


hatched area. 
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Fig. 5. The proton and He nucleus (x5) energy spectrum during 1959. 

The data are from Webber and McDonald [1964], and Freier 
and Waddington [1965], The smooth curve that fits the 1963 
data is shown for comparison. 

Fig. 6. The relative proton to He nucleus ratio inside the solar system 
as a function of the kinetic energy per nucleon according to 
the solar wind modulation theory, and with the choice of the 
distribution of the scale size of the magnetic inhomogeneities 
given in the text. See also the text for the definition of the 
parameters K L and ^ 0 /R q . 

Fig. 7a. The proton and He nucleus (x7.2) rigidity spectra during 1965. 
The data are the same as in Figure 2. 

Fig. 7b. The proton and He nucleus (x7.2) rigidity spectra during 1963. 
The data are the same as in Figure 2. 

Fig. 7c. The proton and He nucleus (x7.2) rigidity spectrum during 1961 
The data are the same as in Figure 2. 

Fig. 7d. The proton and He nucleus (x7.2) rigidity spectrum during 1959 
The data are the same as in Figure 5. 
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